Translocation of a Polymer through a Nanopore across a Viscosity Gradient 
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The translocation of a polymer through a pore in a membrane separating fluids of different 
viscosities is studied via several computational approaches. Starting with the polymer halfway, we 
find that as a viscosity difference across the pore is introduced, translocation will predominately 
occur towards one side of the membrane. These results suggest an intrinsic pumping mechanism 
for translocation across cell walls which could arise whenever the fluid across the membrane is 
inhomogeneous. Somewhat surprisingly, the sign of the preferred direction of translocation is found 
to be strongly dependent on the simulation algorithm: for Langevin Dynamics (LD) simulations, 
a bias towards the low viscosity side is found while for Brownian Dynamics (BD), a bias towards 
the high viscosity is found. Examining the translocation dynamics in detail across a wide range of 
viscosity gradients and developing a simple force model to estimate the magnitude of the bias, the LD 
results are demonstrated to be more physically realistic. The LD results are also compared to those 
generated from a simple, one dimensional random walk model of translocation to investigate the role 
of the internal degrees of freedom of the polymer and the entropic barrier. To conclude, the scaling of 
the results across different polymer lengths demonstrates the saturation of the preferential direction 
with polymer length and the non-trivial location of the maximum in the exponent corresponding to 
the scaling of the translocation time with polymer length. 

PACS numbers: 87.15.ap, 82.35Lr, 82.35.Pq 

Keywords: Nanopore, translocation, polymer, MD simulations, Langevin Dynamics, Brownian Dynamics, 
viscosity gradient 



INTRODUCTION 



The transport of polymeric molecules from one space 
to another through a constricting passage represents a 
fundamental process with a diverse set of applications. 
Within biological systems, there are numerous examples 
that fall under this umbrella including the transport of 
DNA, RNA, and proteins across cell membranes the 
packing of DNA or RNA into - and subsequent release 
from - viral capsids [2] , and the passage of proteins across 
the eukaryotic endoplasmic reticulum [3|. In addition to 
these naturally occurring examples, there is considerable 
interest in this topic due to the implications for the de- 
sign of nanofluidic devices. By far the most prominent 
of such applications concerns using the translocation of 
DNA through a synthetic nanopore as the basis for rapid 
and cheap sequencing technologies 

As a consequence, there has been a large volume of re- 
search on this topic including many theoretical and sim- 
ulation studies (c.f. [f| and references therein). In the 
majority of these studies, the conditions across the pore 
are taken to be equivalent with the exception that, for 
the case of driven translocation, there is a driving force 
towards one side. That is, modelling parameters such 
as the quality of the solvent and the membrane-polymer 
interactions are the same on both sides. However, for 
both biological and synthetic applications, these condi- 
tions are unlikely to be uniform across the membrane. 
For example, considering translocation across a cell wall, 
the type and concentration of solutes is not equivalent for 
the intracellular and extracellular fluid and, as an associ- 



ated effect, the viscosity of the fluid is different across the 
membrane. While the majority of the literature neglects 
such effects and considers translocation driven by an ex- 
ternal field, there have been a number of studies consid- 
ering a bias which results from an asymmetry across the 
membrane. Examples include a bias arising from varying 
solvent conditions across the pore 0, 0| , differing concen- 
trations of obstacles across the membrane Q , chaperone 
assisted translocation 0, [l(| and translocation induced 
by the adsorption of monomers to one side of the mem- 
brane [iIIElI. 

In this work, we study translocation in the presence of 
a viscosity gradient where the viscosity on the cis side is 
different than the viscosity on the trans side. Performing 
both Langevin Dynamics (LD) and Brownian Dynamics 
(BD), we find that a viscosity gradient introduces a bias 
to monomers at the interface. Starting with the polymer 
halfway through the pore, this bias establishes a preferred 
direction for translocation which grows with an increas- 
ing discrepancy between viscosities. In what may be a 
surprising result, which side is preferred depends on the 
simulation algorithm: in LD, more events occur to the 
low viscosity side; in BD more events occur to the high 
viscosity side. To characterize this drastic difference of 
results, simulations exploring the details of the dynamics 
are performed. A simple force model is also employed to 
estimate the magnitude of the viscosity gradient bias for 
both LD and BD results. From this analysis, we take the 
LD results to be more physical and additional simulations 
at different polymer lengths N are performed to exam- 
ine the dependence of the scaling of both the preferential 
direction and translocation time r on the magnitude of 
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the viscosity gradient. We find that the strength of the 
preferential direction slowly grows with TV. Comparing 
these results to those obtained for a single random walker 
at a viscosity interface indicates that this N dependence 
arises from the internal degrees of freedom of the poly- 
mer. For the translocation time, the scaling exponent a 
obtained from r ~ N a is found to decrease as the bias at 
the interface increases, as expected. However, the max- 
imum in a is shown to occur not when the viscosities 
across the pore are equal, but rather when the viscosity 
on one side is slightly higher than the other. The compet- 
ing mechanisms which lead to this result are discussed. 

SIMULATION SETUP 

To model the polymer, we employ an approach that is 
by now common for coarse-grained simulations 13[. To 



prevent overlap of the monomers, a shifted and truncated 
Lennard- Jones potential (often called the WCA potential 
[3]) is used for excluded volume interactions. Defining 
e to be the energy of the interaction and a to be the 
monomer diameter, the potential as a function of the 
centre to centre distance between monomers, r, is given 
by 



Uwca{t) 




e for r < r c 
for r > r r 



(1) 



where r c = 2 1 / 6 cr is the cut-off distance. Bonds between 
monomers along the polymer are modelled by the FENE 
potential which is given by 
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(2) 



Following the model of Kremer and Grest [15| we set 
k = 30e/cr 2 and ro = 1.5a. A continuous surface is used 
to model the nanopore-containing membrane. To achieve 
a pore that is as large as possible while still ensuring a 
sing le file process, the radius of the pore is set to 1.5er 

m 

In this work, an explicit solvent is not included and 
instead the effects of the solvent are included implicitly 
in the equation of motion for the monomers. While this 
approach neglects hydrodynamic interactions, the com- 
putational savings are significant and allow for a much 
more thorough investigation. Two simulation approaches 
with implicit solvent are employed: Langevin Dynamics 
(LD) and Brownian Dynamics (BD). In LD, the effects 
of a solvent are included implicitly by adding two terms 
to the equation of motion. First, a drag force Ff propor- 
tional to the monomer velocity v represents the dissipa- 
tion aspects of monomer- solvent interactions. Second, a 
random force models the fluctuations. The final equation 
of motion is given by: 



where m is the mass of the monomer, U(r) is sum of the 
conservative potentials, ( is the friction coefficient, and 
R(t) is the random term. To ensure that the fluctuation- 
dissipation theorem is obeyed, R(t) satisfies the following 
criteria: 

(Ri(t)) = (4) 
(R\(0) ■ Rjit)) = 2k B TQ§{t)5 ij (5) 

where the subscript i,j denote components along carte- 
sian coordinates. 

Considering the dissipation term, the drag on the par- 
ticle moving through a fluid is given by Stoke's relation 



-6irr/R p v, 



(6) 



where R p is the radius of the particle. The friction coef- 
ficient is thus proportional to the viscosity, £ = 67rr)R p , 
and considering a fixed particle size, varying £ thus di- 
rectly corresponds to varying the viscosity of the fluid. 
To make this identification, £ has units of \J em /a 2 and 
we thus define a dimensionless viscosity given by 

*=-r^Ti' (7) 

In BD, the dynamics are taken to be in the overdamped 
limit where the inertia of a particle is negligible com- 
pared to the drag and random forces. A particle thus 
immediately reaches terminal velocity and the trajectory 
consists of a series of uncorrelated jumps induced by the 
random force. In these simulations, Eqn. [3] is modified 
by dropping the inertial term to yield: 



(v = -VU(r)-R(t). 



(8) 



mf = -VU(r) -(v + R(t), 



(3) 



Performing both LD and BD, the same model for the 
polymer and nanopore is used between simulations. How- 
ever, in BD, the WCA potential is capped at 75 e/a to 
prevent breaking of the bonds due to large overlaps of 
monomers arising from large jumps. 

The system is constructed by having different viscosi- 
ties of fluid on either side of the membrane. As shown 
in Fig. [TJ the interface is placed halfway through the 
pore. Particles crossing the surface are treated as point 
particles such that the friction coefficient in eq. [3] is cho- 
sen based on which side of the interface the centre of 
the monomer is on. As we study unbiased translocation, 
simulations begin with the polymer halfway through the 
pore with an equal number of monomers on cis and trans. 
The middle monomer, which is considered to be on cis 
at t = 0, is fixed to allow the polymer to equilibrate. 
Following equilibration, the polymer is released and the 
direction and time of translocation are recorded. For 
all simulations, the viscosity on the cis side is held at 
i)c = 1.0. The viscosity on trans is varied from t)t = 0.1 
to 9.0. 

For the LD simulations, polymers of length N = 
25,49,75,99 are studied. In the BD simulations, a poly- 
mer of length N = 49 was simulated. 
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FIG. 1: Schematic of unbiased translocation across a viscosity 
gradient. The viscosity on the cis side is held fixed at 77c = 
1.0 while the viscosity on the trans side is varied from ?7t=0.1 
to ?jt=10. 



RESULTS 



so 



70 



in 
c 

9> 60 

CD 

so 
-t— • 

C 

1 40 

CD 



20 



1 jS**^ 


-^events to cis 








- 


(F v. 
J 'S. 

P V. 




-/ »» 

— 


- 

--^events to trans 


1 1 1 1 1 


1 1 1 



1 2 3 4 „ ? 6 7 S ! 

FIG. 2: Fraction of events to the cis side (solid) and trans 
side (dashed) as a function of fjr for LD simulations with 
a polymer of size N=^9. The low viscosity side is always 
preferred. 



Langevin Dynamics 

Beginning with the results from LD simulations, the 
data for the preferential direction is displayed in Fig. [5] 
For f]T < f)Ci more events occur towards trans. Con- 
versely, for t)t > rjc, more events occur towards cis. 
Hence, the low viscosity side is always preferred (with 
equal probability at fjx — fjc as required. Further the 
strength of the preferential direction increases with an 
increasing viscosity gradient; for fjr/fjc — 9) nearly 80% 
of all events end up on the cis wall. 

The respective translocation times are shown in Fig. 
[31 Unsurprisingly, the translocation time to the low vis- 
cosity is always less than the time to the high viscos- 
ity side. However, the difference between tt and tc is 
small compared to the variation of both with ryr- That 
is, although increasing t)t introduces a preferred direc- 
tion towards cis, the net time increases due to the slower 
dynamics on the trans side. 



Brownian Dynamics 

The results for the preferential direction obtained via 
BD simulations are shown in Fig. In exact disagree- 
ment with the LD results, the high viscosity result is 
always preferred in BD. Further, the strength of the pref- 
erential direction grows much faster with t)t than in the 
LD results. Here, when fjr/fjc — 2, essentially all events 
end up on the trans wall. 

The corresponding translocation times are shown in 




FIG. 3: Translocation time to the cis side (solid) and trans 
side (dashed) as a function of f\T for LD simulations for a 
polymer of size iV=49. The inset highlights the behaviour at 
low viscosities. 



Fig. [5] Again, significant disagreement with the LD 
results is found. In LD, both tt and tq increased mono- 
tonically with increasing ffr. In the BD results, tt and 
tc are maximum in the vicinity of fj-r/fjc = 1 and the 
translocation times on either side decrease with an in- 
creasing viscosity difference. 
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FIG. 4: Percent of events to the cis side (solid) and trans side 
(dashed) as a function of fjr for BD simulations for polymers 
of size JV=49. The high viscosity side is always preferred. 
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FIG. 5: Translocation time to the cis side (solid) and trans 
side (dashed) as a function of fjr for BD simulations at a 
polymer size of iV=49. Data is shown only when at least 10 
events are recorded. 



LD vs. BD 

Diffusion in an inhomogeneous medium 

It is worth discussing the marked difference between 
the LD and BD results. We first note that both results 
follow from work we have recently published studying 



a single random walker at a viscosity interface 17j. In 
that work, it was shown that for BD simulations, par- 
ticles near the interface will tend to accumulate on the 
high viscosity side while for LD simulations, the particles 
preferentially end up on the low viscosity side. 

The difference fundamentally amounts to the extent to 
which the particle "feels" the interface as it is crossing it. 
In BD, there is no memory in the system; the dynamics 
at each time step are independent of the dynamics at all 
previous time steps. Correspondingly, a particle that is 
in the low viscosity side will jump over the interface - 
with a jump length given by the low viscosity - and land 
in the high viscosity side. Thus, a particle going from 
low to high viscosity will jump far into the high viscosity 
side. Conversely, a particle jumping from high to low 
viscosity regions will do so with a short jump length and 
land near to the interface. The net effect is that there is 
a bias at the interface favouring particles moving to the 
high viscosity side. Given the independent dynamics, the 
particle does not see the interface as it is crossing it - it 
simply lands in a region of different viscosity. 

In LD simulations however, there is memory in the sys- 
tem due to the inertial term in the equation of motion. 
That is, there is a finite correlation time in the system 
and it takes a measurable amount of time for the veloc- 
ity of the particle at any given time to be washed out. 
Correspondingly, as a particle is crossing the interface, it 
can feel the interface since the rate at which the particle 
inertia is being damped will change. This means that 
particles crossing from low to high will be stopped short 
(compared to dynamics in the low side) while particles 
crossing from high to low will land further past the in- 
terface (compared to dynamics in the high side). This 
difference amounts to an effective bias at the interface 
favouring the low viscosity side. Hence, in BD, particles 
tend to accumulate on the high viscosity side while in 
LD, there is a bias favouring the low viscosity side [l7j |. 
To be more precise, BD corresponds to an Ito formu- 
lation and LD corresponds to an isothermal formulation 
for the problem of diffusion in an inhomogeneous medium 

Applying these ideas to polymer translocation gives 
context to the contradictory results shown above. In LD, 
monomers at the interface will tend towards the low vis- 
cosity side and, correspondingly, a net bias to the low 
viscosity side is established yielding the preferential di- 
rections shown in[2j For the BD results, particles tend to 
accumulate on the high side and translocation to the high 
side is the most likely outcome. A comparison between 
the BD and LD dynamics of a particle at the interface is 
shown in Fig. [6] In both cases, we consider a monomer 
at the interface jumping in the preferred direction (high 
viscosity for BD, low viscosity for LD). The jump to the 
right or left yields an increase in chain tension with re- 
spect to the monomers left behind. In BD, the particle 
moves slower once it is on the high viscosity side and it is 
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effectively trapped. To resolve the increased tension, the 
chain on the low viscosity side - which moves and relaxes 
more freely - will be pulled towards the high viscosity 
side. Hence, for BD, both the bias and resulting chain 
tension favour translocation to the high viscosity side and 
this results in the rapid establishment of a preferential di- 
rection as seen in Fig. |4j In LD, the particle jumps to the 
low viscosity side. However, as the monomers remaining 
in the high viscosity move slowly, there will be increased 
tension tending to bring the monomer back to the inter- 
face. For LD, the bias and resulting chain tension thus 
work against each other and, as seen in Fig. [2j the estab- 
lishment of a preferential direction is much more gradual 
than in BD. 
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FIG. 6: Schematic for the dynamics of monomers crossing the 
viscosity interface for BD and LD. 



Measuring the dynamics 

To get insight into the different behaviour for the 
translocation times, we use the physical picture outlined 
in Fig. IS] to develop a simple test. Roughly speaking, 
there are two factors which affect the translocation time: 
the amount of diffusive motion and the time scale of the 
dynamics. The establishment of a bias at the interface 
would be expected to reduce the amount of diffusive mo- 
tion and thus speed up translocation. Conversely, if the 
bias is generated by increasing the viscosity on one side 
of the membrane, the reduced rate of dynamics would be 
expected to increase the translocation time. To investi- 
gate these opposing contributions for both BD and LD, 
simulations were performed were the translocation coor- 
dinate s was recorded at each time step. Here, s is defined 
to be the number of monomers on trans or, equivalently, 
the monomer currently in the pore where the numera- 
tion begins with s — at the end of the polymer on the 



trans side. Taking a specific monomer, denoted s*, two 
tests were performed. First, the number of times that 
the polymer is displaced by a distance As was recorded 
by measuring the number of times that monomers fur- 
ther along the chain s* + As and back along the chain 
s* — As were in the pore. Second, the time required to 
achieve the displacement was recorded. The first test 
thus indicates the amount of diffusive motion as it indi- 
cates the number of times the polymer travels back and 
forth through the pore. The second test measures the 
time scale of the dynamics. 

Note that as noise is directly added to the equation 
of motion for each monomer, there will always be sig- 
nificant local fluctuations. As we are interested in the 
dynamics on a large time/length scale, a slightly larger 
displacement is chosen and we set As = 3. The polymer 
length is set to 49 and we concentrate on the results for 
Vt > fjc- Further, only trajectories ending at the pre- 
ferred side are considered (in fact, as much of this work 
will be done at a 100% preferential direction, this is a 
necessity since events occur only in one direction). For 
BD, this means events occurring to trans while for LD 
this means events to cis. In either case, the displace- 
ments in the direction of the preferential direction are 
considered positive while those in the other direction are 
negative. For similar reasons, calculations are performed 
only for the monomers initially on the opposite side of 
the membrane as it is necessary for these monomers to 
pass through the pore in order for translocation to be 
achieved. Hence, for BD, the number of events forward 
and backward are averaged over s* — 25 to s* = 46; for 
LD over s*=23tos*=2. For each s*, the number of 
times that s = s* +3 is recorded as a N + for BD (and N_ 
for LD); the number of times that s = s* — 3 is recorded 
as a N- for BD (and N + for LD). Associated times for 
these events are given by t + and f _ . Results for the BD 
simulations are shown in Fig. [7J 

In Fig. [7] a, note the rapid decrease in N + and iV_ 
as t\t increases. By tjt — 10, the values have saturated 
at N + = 1 and iV_ = 0. That is, the polymer never 
goes backwards (at least as far as 3 monomers) and goes 
forwards only once: there is no diffusive motion at the 
interface. While the timescale for the forward motion t+ 
is increasing with t]t as expected, the reduction of iV_|_ 
is dominant and, as shown in Fig. [7]b, the translocation 
time plummets to less than 1 /4 of its initial value by r/x — 
10. In fact, by t\t — 10, tt has started to increase as the 
reductions in N + , iV_ have saturated while t + continues 
to increase. 

The corresponding results for LD are shown in Fig. [5] 
In Fig. [8] a, contrary to the BD results, N + and N- ini- 
tially increase. A physical picture for this result is given 
in Fig. [6l Although monomers tend towards the low vis- 
cosity side, they are pulled back from the slow moving 
chain on the high viscosity side yielding increased fluctu- 
ations. There is agreement with the BD results in that t + 
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FIG. 7: Dynamics of monomers at the interface for BD for a 
polymer of size N = 49. a) displays the number of displace- 
ments forward N+ and backwards N- along with the average 
times for these events given by t+ and t-. b) displays the 
average net trans translocation time. 



FIG. 8: Dynamics of monomers at the interface for LD for a 
polymer of size N = 49. a) displays the number of displace- 
ments forward 7V+ and backwards TV- along with the average 
times for these events given by t+ and t_. b) displays the 
average net translocation time and the number of events oc- 
curring towards the low viscosity side. 



increases with increasing t\t as expected. Consequently, 
as both N + and t + increase, the net translocation time 
initially increases as shown in Fig. [8] b. The initial in- 
crease in N + is short lived and, after rjx ~ 5, N + decays 
with increasing rjx in agreement with the BD results. 
However, contrary to the BD results, the decay of N + 
is weaker than the increase in t+ and the translocation 
time continues to increase. At even higher viscosities, the 
balance of these factors reverses and after t]t = 30, tq 
slowly decreases with increasing ryr- As shown in Fig. [8] 
b, the peak in tc coincides with the preferential direction 
plateauing at 100%. 



Estimating the bias due to the viscosity gradient 

To measure the bias at the viscosity interface for both 
the LD and BD simulations, a simple force model is em- 
ployed. In this approach, the viscosity is held uniform 
a t t)c = f)T = 1 and an external force is applied to the 
monomers at the interface by applying a force F to any 
monomer which is in the pore. We again focus on the 
/V=49 case and start with the polymer halfway. Perform- 



ing these simulations for values of F ranging from 0.01 to 
10, the preferential direction and translocation time are 
recorded. Figure [§] displays the results for the preferen- 
tial direction results for the force model applied to both 
LD and BD simulations along with the results for the 
BD viscosity gradient simulations. In Fig. |H1 there is ex- 
cellent agreement between the external force simulations 
from LD and BD approaches. This result demonstrates 
that, for a system of uniform viscosity, LD and BD simu- 
lations generally do give equivalent results (although not 
shown, good agreement between the translocation times 
is also found with only a uniform rescaling of the times 
between the approaches being required). Good agree- 
ment between the viscosity gradient case and the exter- 
nal force cases is obtained by rescaling the force such 
that Afj = t]t — t]c — 1-5-F. The preferential direction 
results thus allow us to estimate the magnitude of the 
bias at the interface induced by the viscosity gradient via 
F = Afj/ 1.5. Considering that typical values for external 
fields in simulations examining driven translocation are 
in the range of 0.5 - 3 (in equivalent units), this result 
indicates that the bias due to the viscosity gradient in 
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FIG. 9: Matching of the preferential direction for the BD 
viscosity gradient results (red, events to the high viscosity 
side), the external force model applied to BD simulations at 
a uniform viscosity (blue, events in the direction of the force), 
and the external force model applied to LD simulations at a 
uniform viscosity (green, events in the direction of the force). 
The polymer size is set to iV=49. Good agreement is found 
taking Afj = 1.5F. 



BD simulations is of considerable magnitude. 

The same data plotted along with the LD viscosity gra- 
dient results is shown in Fig. QUI Here, to get agreement 
between the external force and viscosity gradient results, 
the force must be rescaled by Af] — 135F indicating that 
F = Ary/135. The force induced by the viscosity gradi- 
ent in LD simulations is much smaller than that typically 
employed to examine biased translocation. Similarly, it 
is nearly two order of magnitude smaller from that gen- 
erated in BD simulations. This discrepancy also demon- 
strates how much more sensitive the BD simulations are 
to viscosity changes with a much stronger preferential di- 
rection being established at smaller values of Afj. Note 
that the value of 135 seems to be independent of the 
polymer length N. For simulations performed at shorter 
polymers (N — 25) and longer polymers (N = 75), using 
a factor of 135 was found to give good agreement between 
the viscosity gradient and force model results (data not 
shown) . 

The inset to Fig. [TU] displays the results from "force 
balance" simulations in which the force is applied in the 
opposite direction to the effective force arising from the 
viscosity gradient. Relating the external force to the vis- 
cosity gradient by A77 = 135F, the forces effectively can- 
cel each other out and the preferential direction all but 
vanishes. This result further corroborates the external 



FIG. 10: Matching of the preferential direction for the LD vis- 
cosity gradient results (red, events to the low viscosity side), 
the external force model applied to BD simulations at a uni- 
form viscosity (blue, events in the direction of the force), and 
the external force model applied to LD simulations at a uni- 
form viscosity (green, events in the direction of the force). 
The polymer size is set to iV=49. Good agreement is found 
taking Afj = 135F. The inset shows the preferential direction 
when the external force and the bias due to viscosity gradient 
oppose each other. Using the Arj = 135F relation, the pref- 
erential direction is effectively nullified across the viscosity 
range studied. 



force model as well as the estimate of the magnitude of 
the bias due to the viscosity gradient. 

Matching the preferential directions allows us to es- 
timate the magnitude of the bias at the interface due 
to the viscosity difference. In comparing the transloca- 
tion times, there is another factor we must consider. For 
the external force cases, an increasing F necessarily de- 
creases r. However, in the case of the viscosity gradient, 
as discussed earlier, there are two factors to consider: i) 
the reduction in r due to an increased bias and ii) an 
increase in r as the average viscosity in the system is be- 
ing increased. To compare r between the external force 
and viscosity gradient simulations, we must account for 
this latter factor. To do so, r from the viscosity gradient 
simulations must be rescaled by a factor that depends on 
the viscosities across the membrane: g(f)c,fjr)- We have 
recently shown that although the translocation time is in- 
dependent of viscosity at very low fj values, it decreases 
as I/77 for intermediate to high viscosities (i.e., fj > 1). 
Using this result and recalling that in these simulations 
we fix fjc and vary f/T, we assume the following form for 
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the g(fjc,Vr) scaling factor: 

givcVr) = Mm - fjc) + f)c- (9) 

This form captures the dependence on the viscosity dif- 
ference Afj = fjr— fjc and ensures that, in the limit where 
Afj = 0, the inversely proportional scaling of the translo- 
cation time with uniform viscosity given by g = fjc — fjr 
is obtained. 

Beginning with the BD case, Fig. QT] displays the 
translocation time corresponding to the external field 
model te plotted against 1.5F. Here we show the results 
from the driving force applied to the BD simulations such 
that t e {F = 0) = T Af j(Arj = 0). Fig. [TTIalso shows the 
translocation time corresponding to the viscosity gradi- 
ent results ta^ where the time has been normalized by 
g. Good agreement is found when A = 1 such that 

g = (fjr - fjc) + fjc (10) 
= VT. (11) 

This indicates that in BD where the events occur to the 
high viscosity side, the rescaling of the translocation time 
to account for the changing viscosity in the system is 
achieved simply by normalizing r by the viscosity on the 
high viscosity side. This simple result demonstrates that 
as events occur towards this side, motion of the entire 
polymer through the high viscosity side is the rate lim- 
iting step for translocation. Hence, for the BD results, 
good agreement for both the preferential direction and 
the translocation time is achieved by matching the ex- 
ternal field to the viscosity gradient via Afj = 1.5F and 
normalizing the time scale of the dynamics by f\T- 

Figure [TT] also demonstrates that there is a transition 
from an essentially F independent region to a region 
where r decreases as 1/F around F* = 0.1/1.5 = 0.667. 
This crossover reflects the transition between the dynam- 
ics being determined by diffusive or driving processes 
and has been found before for studies examining driven 
translocation at low forces [2!|. Below F*, diffusion is 
dominant and, since the external force plays a minor role 
in determining the translocation time, r is relatively in- 
dependent of F. Above this point, the driving force dom- 
inates and r decreases proportional to F. For the viscos- 
ity gradients studied, the BD results correspond to the 
high field limit again reflecting the large impact of the 
viscosity gradient in BD simulations. 

For the LD case, Tp plotted against 135F and r v /g 
plotted against lS.fi is shown in Fig. [12] In this plot, 
the results for the external field generated from LD sim- 
ulations are shown. In the LD case, good agreement is 
found when A = 0.34 such that 

g = 0M(fj T -fjc)+fjc (12) 
= 0.65(j? r - 1) + 1. (13) 

Contrary to the BD results where a direct scaling of 1 /fjr 
was obtained, in LD there is a weaker dependence of the 
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FIG. 11: Matching of the translocation times by plotting the 
BD viscosity gradient times normalized by r\T along with the 
BD external field times. The polymer size is N = 49. 

translocation time on the magnitude of the viscosity on 
the high viscosity side. This result reflects that in LD, 
the events occur towards the low viscosity side and thus 
the impact of increasing the viscosity on the other side is 
reduced compared to BD. 
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log(Afj),log(135F) 

FIG. 12: Matching of the translocation times by plotting the 
LD viscosity gradient times normalized by 0.65(?7t — 1) + 1 
along with the LD external field times. The polymer size is 
TV = 49. 

Figure IT21 also shows that, contrary to the BD results, 
some of the LD results lie in the region where r is rel- 
atively unaffected by F. Again, the viscosity gradient 
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results from LD correspond to low fields while from BD 
they correspond to high fields (this is particularly true re- 
calling that the range of Afj simulated for BD was smaller 
than that for LD). 

Hence, for both BD and LD, a force equivalent to the 
bias introduced by the viscosity gradient can be obtained. 
For BD, the magnitude of the force is nearly 100 times 
larger than that for LD. In fact, as the equivalent BD 
force is so large, the reduction in the translocation due 
to having a preferential direction is much greater than the 
increase in the translocation time due to having a higher 
average viscosity in the system. Hence, the translocation 
time decreases with an increasing t)t. The converse is 
true in LD where the equivalent force is quite weak and 
the reduction in r due to the preferential direction is 
negligible compared to the increase in the time scale of 
the dynamics due to increase fjr- Correspondingly, r 
increases as t)t does. 



Physicality of the alogrithms 

The drastic differences between these results raise the 
question as to which one is physical. We believe the 
LD result is the correct one for this scenario. First, the 
equation of motion for LD is more correct in including 
the inertial term since it does not require the assumption 
of overdamped dynamics implicit to BD. Further, the 
results for the translocation time for the BD simulations 
seem to be unphysical. It is difficult to believe that, given 
an average translocation time when t)t — fjc = L the 
translocation time is reduced by a factor of four when the 
viscosity on the trans side is doubled. As demonstrated 
by Fig. the BD approach quickly yields a ratchet in 
which monomers jump from low to high viscosity and 
never cross back. In fact, these dynamics result from the 
monomer jumping over the interface without feeling it 
such that a monomer from the low viscosity side is essen- 
tially artificially "injected" into the high viscosity side 
without experiencing the change in viscosity. As shown, 
this quickly suppresses diffusive motion at the interface, 
yields a rapid establishment of a 100% preferential direc- 
tion, and produces a dramatic drop in the translocation 
time. These results are consistent with having a pump 
for monomers from the low viscosity side to the high vis- 
cosity side. As there is no source for the energy required 
to drive such a pump, the BD results seem unphysical 
indicating that at least this implementation of BD is in- 
appropriate for the system under study. 

As the LD simulation results are taking to be correct 
for this scenario, the remainder of the paper focuses on 
the LD results. However, it is worthwhile to bring at- 
tention to the remarkably different results obtained be- 
tween LD and BD as this example demonstrates that 
while these approaches generally yield the same result, 
scenarios exist where the difference between these ap- 



proaches can dictate the final result. 



One dimensional random walker in a viscosity 
gradient 

As mentioned, the LD results follow from a study we 
have published concerning a single particle at a viscos- 



ity interface 17[. In that work, the particle was initially 



placed at a sharp interface equidistant between two walls. 
Here, we build on this approach to make it more analo- 
gous to the case of polymer translocation. Comparison 
between the simplified model and the polymer results 
will then give insight into the role of the internal degrees 
of freedom of the polymer and the effect of the entropic 
barrier. In the aforementioned study, Monte Carlo ap- 
proaches for simulating the system were introduced and 
tested along with Langevin and Brownian Dynamics sim- 
ulations. While these approaches can be extended to the 
case of polymer translocation studied here, in the current 
article we present results only from Langevin dynamics 
simulations for the sake of simplicity of comparison to 
the polymer results. 

To study translocation, several studies (including the 
earliest theoretical treatments of the problem) have used 
a quasistatic approximation to reduce the process to a 
one dimensional random walker (1DRW) IB-HI- The 
basic idea is that, if the polymer can be considered in 
equilibrium at all times, translocation can be represented 
as a single random walker traversing an entropic bar- 
rier which arises from confining the polymer within the 
nanopore. In a similar fashion, the effects of the viscos- 
ity gradient can be included in this model. Considering 
a polymer halfway through the pore, there will be N /2 
monomers on cis experiencing i)c and N/2 monomers on 
trans experiencing fjT- The average friction coefficient for 
this configuration is thus (?7c+?7t)/2. Defining a translo- 
cation coordinate s which indicates how many monomers 
are on the trans side, this definition can be extended: 



rj(s) 



(N - s)fj C + sr) T 
N 



(14) 



where N , being the polymer length, corresponds to the 
separation between the absorbing walls for the 1DRW. 

Using this form, Langevin dynamics simulations for a 
one dimensional random walker in a linear viscosity gra- 
dient were performed. Simulations were done for both 
a "free" particle (subject only to the viscosity gradient) 
and a particle traversing an entropic barrier. The form 
for the entropic barrier is given by the quasistatic ap- 
proximation 2ol- 25j| : 



U 



( S 



-kTln 



1 



(15) 



where 7 is the surface exponent set to 0.69 here [26[ such 
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that a force given by 

lisl = 1 

kT N 
is applied to the particle. 



1 — — 

1 N 



N 



(16) 




FIG. 13: Schematic of ID random walker in a viscosity gra- 
dient 

As shown in Fig. [T3l the particle was initially placed 
in between two absorbing walls. Simulations were per- 
formed for fjx = 0.1 to 9.0 with f)c held at 1 where fjr 
is now the viscosity at the right wall and fjc is the vis- 
cosity at the left wall. The percent of events absorbed 
at either wall and the mean first passage time (MFPT) 
were recorded. The results for the preferential direction 
are shown in Fig. [TJ] along with the polymer results. 
In agreement with the LD polymer results, the particle 
preferentially ends up at the low viscosity wall. 

Comparison between the 1DRW and polymer results 
indicate several things. First, the consistent establish- 
ment of a preferential direction to the low side for all 
three simulation approaches demonstrates that the in- 
ternal degrees of freedom of the polymer are not needed 
to obtain this result; a preferential direction can be es- 
tablished for just a single particle with or without the 
cntropic barrier (this is also true for the case of a sharp 
interface 

However, the internal degrees of freedom of the poly- 
mer do have an affect on the dynamics. Examination 
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FIG. 14: Percent of events absorbed at the cis wall (solid) and 
the trans wall (dashed) for the polymer simulations (red), the 
1DRW viscosity gradient simulations (blue), and the 1DRW 
simulations including the entropic barrier (green). The poly- 
mer size is set to N—49. 



of the preferential direction plotted logarithmically in 
fj shows that the results for the 1DRW are symmetric: 
n-T /ncifir) = tic/titO-/Vt)- Conversely, the results for 
the polymer translocation are clearly not symmetric. The 
difference between these two cases is that for the poly- 
mer, the internal degrees of freedom introduce another 
time scale into the system: the time required for the poly- 
mer to relax. The relaxation time breaks the symmetry 
between fjc/fjr and fjr /fjc- We have recently published 
a study examining the unbiased translocation dynamics 
as a function of viscosity. These results indicate that, for 
N=A9, the dynamics at fj = 0.1 are essentially quasistatic 
while for fj = 1 and higher, nonequilibrium effects would 
be expected. Hence, at fjr /fjc = 0.1, the dynamics on 
trans are quasistatic while those on cis are not. However, 
at ffr/f]c = 10, neither the dynamics on trans or cis are 
expected to be quasistatic. The result is the asymmetric 
results shown in Fig. [T4l 

Interestingly, the nonequilibrium effects appear to be 
strengthening the preferential direction since ncifir) > 
nr(l/fjT)- To give context to this, consider a scenario 
in which the viscosity on both sides was low enough to 
yield quasisatic dynamics. For this case, the preferential 
direction would arise from the bias felt by monomers at 
the interface towards the low viscosity side. If the vis- 
cosity on both sides is now raised (while maintaing the 
same ratio) such that neither side corresponds to qua- 
sistatic dynamics, another bias is introduced due to the 
differing relaxation times across the interface. Given that 
the relaxation time will be lower on the low viscosity 
side, movement of monomers towards the low viscosity 
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side will in general be more successful than attempts to 
move monomers to the high viscosity side (for a detailed 
study of the effects of viscosity on memory effects during 
translocation, see 27] ). This extra bias resulting from 



this discrepancy of relaxation times adds to the tendency 
for monomers towards the low viscosity side and thus 
yields a stronger preferential direction. 

Comparing the 1DRW results, a noticeable drop in the 
strength of the preferential direction is observed when 
the effects of the entropic barrier are included. For this 
case of a continuous viscosity gradient, one can consider 
the particle to be experiencing a bias to the low viscos- 
ity wall at all locations. In the absence of an entropic 
barrier, the net effect of this bias yields a strong prefer- 
ential direction. When the entropic barrier is added, this 
bias still dominates near the middle of the system where 
the entropic barrier is flat and thus the resulting force is 
negligible. However, when the particle is quite near to 
a wall, the entropic barrier yields a large force towards 
that wall. In these regions, the force arising from the 
viscosity gradient will be negligible. The particle thus ef- 
fectively sees a reduced viscosity gradient corresponding 
to the viscosity drop over the middle of the system. As 
shown previously, a reduced viscosity gradient yields a 
weaker preferential direction and thus there is a drop in 
the strength of the preferential direction on including the 
viscosity gradient in the simulations. This result will be- 
come important when considering the results for different 
polymer lengths discussed below. 
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FIG. 15: Mean first passage time to cis (solid lines) and trans 
(dasehd lines) for polymer translocation (red), a free 1DRW 
(blue), and a 1DRW traversing an entropic barrier (green). 
For ease of comparison, all times are normalized to be equal 
at i)c = t)t- 

An additional bias to the low viscosity side for higher 
f)T may also explain the discrepancies for the times shown 



in Fig. [15] When the times are normalized to be equal at 
Vc = Vt, the 1DRW MFPTs overestimate the polymer 
translocation times for t)t > 1. This result is consistent 
with the polymer experiences an additional bias towards 
the low viscosity side such that r does not grow as fast 
with t)t as for the 1DRW simulations. 1 In contrast, 
the 1DRW results underestimate the translocation time 
when f]T < 1 (inset to Fig. II 5|) . In reducing the problem 
to a one dimensional process, many aspects of transloca- 
tion are neglected. In the 1DRW model, lower viscosity 
values directly correspond to higher diffusion coefficients 
and thus r continues to decrease with decreasing r\T- On 
the other hand, for a polymer in three dimensions, no 
matter how low the viscosity, the polymer must thread 
through the constricting nanopore in order to translo- 
cate. In fact, we have recently shown that at very low 
(uniform) viscosities, the translocation time is indepen- 
dent of the viscosity; friction with the pore is entirely 
dominant and lowering the viscosity does not speed up 
the process [27| ■ Friction with the pore can thus limit the 
reduction of the translocation time and as this effect is 
missing in the 1DRW model, lower MFPTs are predicted. 



Scaling with N 

While the above discussion concerns only the results 
for iV=49, simulations were also performed for N — 
25,75,99. The results for the preferential direction are 
shown in Fig. [TB] The strength of the preferential di- 
rection consistently increases with increasing molecular 
weight. This increase is highlighted in Fig. [T7] where 
the preferential direction is plotted against the poly- 
mer length for the two extreme viscosity gradient cases: 
f)T = 9, fjc = 1 and fjr = 0.1, fjc = 1. For the t)t = 9, 
the number of events towards cis increases to near 90% 
at N = 99. Correspondingly, the saturation of the pref- 
erential direction is evident as the curve must eventually 
plateau at 100 %. For the f\T = 0.1 data, the increase in 
the strength of the preferential direction is more modest 
and the saturation is not evident over the polymer length 
range studied. Comparing between the curves, the dis- 
parity between the curves for viscosity ratios which are 
nearly the inverse of each other again demonstrates the 
asymmetry of the dynamics intrinsic to the polymer re- 
sults. 

Considering the increase in the preferential direction 
with TV, we note that the preferential direction for the 
1DRW is independent of the separation of the walls (re- 
sults are not shown). Beyond trivially increasing the 
MFPTs, there is no effect of increasing the separation 
between the walls. Recalling that the internal degrees of 
freedom of the polymer yield an additional bias towards 
the low viscosity side, this strengthening of the preferen- 
tial direction with increasing N can be explained by not- 
ing that this extra bias will increase with N. That is, as 
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the polymer grows, the disparity in relaxation times be- 
tween the cis and texttrans side grows yielding a stronger 
bias to the low viscosity side and thus a stronger prefer- 
ential direction. 
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FIG. 16: Percent of events towards cis (solid) and trans 
(dashed) wall for polymers of size N = 25 (red), N = 49 
(blue), N — 75 (green), and N = 99 (brown) as function of 
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FIG. 17: Saturation of the preferential direction as a function 
of the polymer size N for the two extreme viscosity gradient 
cases: fjr = 9, fjc = 1 (events towards cis) shown in red and 
f]T — 0.1, fjc = 1 (events towards trans) shown in blue. 

There is another effect to consider here. Recall that 
adding the entropic barrier to the 1DRW model weak- 



ened the preferential direction. As the particle/polymer 
nears translocation, entropic effects dominate the bias 
due to the viscosity gradient and a reduced effect of the 
gradient is observed. It has been shown that for poly- 
mer translocation at viscosities where the process is not 
quasistatic, the effects of the entropic barrier diminish 
with length [28J. That is, entropic effects play a bigger 
role at N=25 than N=99. Correspondingly, the dimin- 
ishment of the viscosity gradient bias due to the entropic 
barrier is reduced as N increases and the strength of the 
preferential direction increases with N. 

Hence, both the growing discrepancy between relax- 
ation times on either side of the pore and the postu- 
late that the polymer experiences a large portion of the 
viscosity gradient result in a preferential direction that 
grows with N. While it may be possible to design tests to 
unravel these effects, for now we simply note that there 
are two possible mechanisms - both arising from the in- 
ternal degrees of freedom of the polymer - which could 
explain this result. 

The scaling of the translocation times was also exam- 
ined. Figure [IS] displays the translocation to both the 
trans and cis sides for TV = 25, 49, 75, 99 at 4 different 
values of the viscosity on the trans side. Note that for 
Vt < fjc , t t < t~c and that for t)t > fjc, t~t > tq- For 
Vt = VCi t t ~ i~c as required. 
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FIG. 18: Translocation times to trans (solid lines) and cis 
(dashed lines) for fjc = 1.0 and fjr = 0.1 (red), 1.0 (blue), 5.0 
(green), and 9.0 (brown) as a function of N. 

Assuming a scaling of translocation time with polymer 
length given by 

r-A Q , (17) 

values for the a exponent can be obtained from the slopes 
of the lines in Fig. [T5] These results are shown in Fig. 
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[T9l There is a maximum just beyond r/r — 1.0. On 
either side of this region, the scaling decreases with an 
increasing difference in viscosity. 
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FIG. 19: Scaling exponent a for trajectories occurring to- 
wards trans (red) and cis (blue) as a function of i)t- 

The physical picture here is that, if translocation 
is predominantly a diffusive process, one would expect 
r ~ A^ 2 (in fact, due to nonequilibrium effects, values 
of a greater than 2 are consistently obtained). On the 
other hand, in the limit where the bias is high enough to 
completely suppress diffusion, one would expect r ~ N . 
Thus, the introduction of a bias reduces the a exponent. 
Considering that the viscosity gradient introduces a bias 
that grows with an increasing difference, a decreases both 
as tt/tc — > and as tt/tc —> oo. 

However, the maximum is not at tt/tc = 1 as this 
argument would imply. Rather, it is just to the right 
of this point. As mentioned, we recently studied the ef- 
fect of viscosity on the unbiased translocation dynamics 
[27l [2^|. It was found that the a exponent increased 
significantly with increasing viscosity due to varying im- 
pacts of nonequilibrium effects. This effect of an in- 
creased a at higher fj thus counters the decrease in a 
as tt/tc — > 00. 

From the data, there is a crossover between which ef- 
fect is dominant: for points just beyond tt/tc = 1, 01 
increases due to the higher viscosity on trans; beyond 
tt/tc — 2, a decreases due to the higher difference be- 
tween f]T and fjc- 



CONCLUSION 

In this work, we examine the translocation of a poly- 
mer through a nanopore when the viscosity differs across 



the membrane. To do so, the viscosity on cis is fixed at 
1.0 and simulations are performed for trans values vary- 
ing from 0.1 to 9. Interestingly, the results were found to 
change drastically between simulation algorithms. For 
the establishment of a preferential direction, the poly- 
mer tends to the low viscosity side in LD while in BD, 
more events occur towards the high viscosity side. For 
the translocation time, r varies essentially linearly with 
the trans viscosity in LD. In BD, r decreases with an 
increasing viscosity difference across the pore such that 
r is a maximum when the viscosity is uniform. 

These conflicting results provide an instructive exam- 
ple of the care which must be taken in choosing a simu- 
lation algorithm when studying a particular system. For 
the study of polymer dynamics, LD and BD generally 
produce equivalent results and, for studies of polymer 
translocation, both have been used extensively. In all 
such cases, the viscosity in the system was uniform and 
thus the results were insensitive to the choice of algo- 
rithm. In fact, in this work for the simple force model 
simulations preformed with a uniform viscosity, BD and 
LD were shown to produce equivalent results. However, 
if we now consider a system where the viscosity on trans 
is twice that on cis, the answer of where a polymer that 
starts halfway will end up turns out to depend on the 
simulation approach: in BD the polymer ends up on the 
"thick" side but in LD the polymer ends up on the "thin" 
side. 

To delve into these results, the details of the dynamics 
and magnitude of the effective bias were studied. In BD, 
increasing the viscosity on trans even a relatively small 
amount corresponds to a considerable force at the inter- 
face. Correspondingly, the driving of monomers across 
the interface quickly suppresses all diffusive motion and 
the translocation time decreases rapidly. In LD, the bias 
introduced from a viscosity gradient is nearly two or- 
ders of magnitude smaller than that produced in BD. 
Likewise, the diminishment of diffusive motion is much 
weaker and the translocation time increases with an in- 
creasing viscosity on trans. From this comparison, we 
believe that, in general, the LD results are correct. Fun- 
damentally, the difference arises from the inertial term 
in the equation of motion in LD which allows monomers 
to sense the interface as they pass over it. In BD where 
there is no inertial term, the dynamics occur entirely ac- 
cording to the viscosity at the initial point of the jump. 
Thus, not only is the former more physically complete, 
but the results are more reconcilable with physical intu- 
ition; it is hard to imagine that doubling the viscosity on 
trans will result in all events occurring to this side in 1 /4 
the time. 

Finally, the scaling aspects of these results were ex- 
amined. We find that the strength of the preferential 
direction increases with increasing polymer length. Com- 
parison to a single random walker at a viscosity interface 
indicates that this N dependence arise from the internal 
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degrees of freedom of the polymer. For the scaling of 
the translocation time with N, the a exponent generally 
decreases as the viscosity gradient grows. This trend is 
expected as the process moves from diffusive to driven. 
However, the maximum a does not occur when the vis- 
cosity is uniform. Rather, if the viscosity on trans is 
increased from one to two, a increases. This result is 
consistent with previous observations that, for a system 
of uniform viscosity, a increases with increasing viscosity. 
There is an interplay between these two effects and thus 
a increases for small increases in the viscosity on trans 
but decreases as the viscosity gradient grows to larger 
values. 

The translocation of a polymer across a viscosity gra- 
dient thus presents a wealth of results. Not only does this 
system represent a clear case where BD and LD yield very 
different results, but the results themselves present im- 
plications for naturally occurring cases of translocation. 
Taking the LD results to be more physical, a polymer 
that is halfway between a thin and thick medium will 
experience a bias driving it towards the thin side. The 
strength of this bias increases not only with the viscos- 
ity gradient, but also with the length of the polymer. 
Likewise, the scaling of the translocation time with the 
polymer length is affected with lower values of a being 
found for larger gradients. Although the simulations indi- 
cate that this bias is a relatively small force, these effects 
are likely to arise in both natural and synthetic cases of 
translocation and thus may be of interest both for recon- 
ciling experimental data with theoretical predictions and 
for guiding the design of nanofluidic devices. 
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